Next-generation sequencing of prostate tumors provides independent evidence of xenotropic murine leukemia virus-related gammaretrovirus contamination.

In 2006, xenotropic murine leukemia virus (MLV)-related gammaretrovirus (XMRV) was isolated from prostate cancer tissue (9). However, subsequent studies have yielded conflicting and controversial results, with widespread detection of XMRV suggested to be the result of contamination with mouse DNA (3). While recent work has focused on PCR-based targeted detection of XMRV (10), the advent of next-generation sequencing (NGS) means that it is now possible to interrogate the entire genomes and transcriptomes of human samples for the unique genomic signatures of thousands of viruses (1, 5, 6). With this in mind, we analyzed whole-genome (DNA-Seq) and transcriptome (RNA-Seq) data from 9 human prostate tumors (6 primary and 3 metastatic), 3 prostate tumor-derived murine xenografts, and 1 benign tissue sample from a pelvic lymph node. The xenograft tumors carried significant amounts of host mouse tissue, thereby acting as positive controls. We mapped all nonhuman RNA-Seq reads to a custom database of 3,932 viral genomes and 1,387 microbial genomes downloaded from the National Center for Biotechnology Information RefSeq database (July 2011) (7) and filtered for reads mapping specifically to MLVs, the family that includes XMRV. Murine genomes are hosts to hundreds of retroviruses (8); consistent with this, the xenograft tumor transcriptomes yielded thousands of MLV reads. Interestingly, 2 primary tumors (P3 and P5) also demonstrated significant enrichment of 5 MLVs, especially XMRV, in their transcriptomes, compared to the other patient tumors (P 1.28 10 8 [t test comparing XMRV reads per million]) (Table 1). To search for evidence of viral integration into the tumor genomes of P3 and P5, we mapped the nonhuman DNA-Seq reads to the same viral database. As expected, the xenograft samples showed enrichment of MLV DNA sequences. However, no DNASeq reads from patients P3 and P5 mapped to any of the MLVs (Table 1). Note that our DNA-Seq human genome coverage averaged 4 , while XMRV is expected to integrate on every chromosome within infected genomes (4). Finally, we interrogated RNA-Seq data for 2 mouse-specific segments of the mitochondrial cytB gene (mouse mitochondrial DNA [mtDNA]; NC_005089.1). The xenograft samples were strongly positive, but tumors P3 and P5 also demonstrated enrichment of mouse mtDNA (Table 1). Indeed, the correlation coefficient across all human samples between the numbers of reads per million mapping to XMRV mtDNA and to mouse mtDNA was 0.982 (note that there was no sequence common to the XMRV mtDNA and mouse mtDNA). To our knowledge, this is first analysis of XMRV in NGS data and provides no evidence for XMRV infection in human prostate tumors. The 2 human samples with enrichment of XMRV reads were also enriched with mouse mtDNA, suggesting a contaminant origin for the viral reads in the samples from those patients. Interestingly, the RNAs from these two samples were extracted at similar time points distinct from those of the other tumor samples and the DNA extractions. The recent emergence of large-scale sequencing projects such as the International Cancer Genome Consortium (2) means that it should soon be possible to mine